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The melting of the Abrikosov vortex lattice in a 2D type-II superconductor at high magnetic 
fields is studied analytically within the framework of the phenomenological Ginzburg-Landau theory. 
Tn , It is shown that local phase fluctuations in the superconducting order parameter , associated with 

rr, ' low energies sliding motions of Bragg chains along the principal crystallographic axes of the vortex 

lattice , lead to a weak first order 'melting' transition at a certain temperature Tm , well below 
the mean fleld Tc , where the shear modulus drops abruptly to a nonzero value. The residual shear 
modulus above Tm decreases asymptotically to zero with increasing temperature. Despite the large 
phase fluctuations, the average positions of Bragg chains at fimite temperature correspond to a 
P3 . regular vortex lattice , slightly distorted with respect to the triangular Abrikosov lattice. It is also 

O ' shown that a genuine long range phase coherence exists only at zero temperature; however, below the 

^ , melting point the vortex state is very close to the triangular Abrikosov lattice. A study of the size 

^ ■ dependence of the structure factor at finite temperature indicates the existence of quasi-long range 

order with S (G) ^ N"^ , and 1/2 < a < 1, where superconducting crystallites of correlated Bragg 
^ ' chains grow only along pinning chains. This finding may suggest a very efficient way of generating 

pinning defects in quasi 2D superconductors. Our results for the melting temperature and for the 
Cu ■ entropy jump agree with the state of the art Monte Carlo simulations. 

~' PACS numbers: 74.60.-w , 74.60. Ge , 74.40. -hk , 74.76.-w 
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I. INTRODUCTION 



Many potentially iniportant|-superconductors , such as some of the High T^ cuprates or the organic charge transfer 
salts K— {BEDT — TTF)2X P), are highly anisotropic compounds with nearly two-dimensional electronic structure. 
^ I These compounds are extremely type-II superconductors with very small inplane coherence length. Consequently, the 
(~\J ■ Ginzburg critical region is relatively large and so one expects drastic deviations from the predictions of the mean field 
\^ \ theory for these materials , due to strong thermal fluctuations in the superconducting order parameter. 
^~^ . In the mixed state at very low temperatures amplitude fluctuations are suppressed , but plpftsje fluctuations can lead 

^N ' to the melting of the vortex lattice at a certain magnetic field Hm (T) , Md < H„i < HcMm A soft shear , Gold- 
^=^ , stone mode, which can be described by long wavelength phase fluctuationaJ is responsible for this remarkable melting 
phenomenon. Unfortunately rigorous analytical approaches to this problem have encountered fundamental difficul- 
ties: laraa-|0rder high temperature perturbation expansion with Borel-Pade approxinrants to the low tenrperature 
j^ ' behaviorQ,El has no indication of an ordered vortex lattice even at zero temperature. The existing non-perturbative 
approaches, have not completely clarified the situation: Renormalization group studieaD , as well as Moirte Carlo 



simulatiorJa have predicted no crystal vortex sta.tp_in a pure 2D superconductor (SC) at finite temperature, while the 

novel functional integral formalism suggested iiilj,E£l, has led to some kind of a vortex liauidJreezing transition without 

C . breaking the U{1) symmetry. Several Monte Carlo simulations, have recently showntJ^tJ that in a 2D , SC a tme 



^ ' vortex lattice melting phase transition takes place at finite temperature and that the transition is of the first orderll3. 
In this paper we present a simple model of the vortex lattice melting in 2D extremely type-II superconductors. 
Our model is based on the observation that at low temperature the main correction to the mean field free energy 
k>( ■ arises from " Bragg chain fluctuations" , namely fluctuations which preserve long range periodic order along a principal 
j_] 1 crystallographic axis in the vortex lattice. This simplification reduces our 2D problem to a ID one , which can then 
d ' be solved exactly. 

Our calculations show that the Abrikosov triangular lattice is subject to strong phase fluctuations. As a result of 
these fluctuations the sharp mean field transition into the Abrikosov lattice state becomes a smooth crossover. Since 
the strength of the phase dependent terms in the SC free energy is relatively small (~ 2% of the SC condensation 
energy) the scale of the crossover temperature, Tcm , is well below the mean field Tc- At temperatures higher than 
Tcm the vortex lattice transforms to an ensemble of strongly uncorrelated vortices, fluctuating independently around 
equilibrium lattice positions. It is found that because of a discontinuous (rotational) symmetry change in the mean 
positions of vortices there is a weak first order transition , superimposed on the smooth solid- liquid crossover, which is 
reflected in a small jump of the vortex system entropy at a certain melting temperature Tm ~ Tcm- Calculation of the 



structure factorj-^lipiiaB that exact long range translational order exists only at zero temperature, in agreement with 
previous resultaj|l3,cll. However, below the melting temperature T^, the vortex state is very close to the triangular 
Abrikosov lattice. Our results for various thermodynamic parameters agree well with numerical calculations in Ref.llil. 

II. THE SLIDING 'BRAGG-CHAINS' MODEL 

Our starting point is the Ginzburg-Landau ( GL ) free energy functional 



Fgl = / d\ 



-«i^(^)i' + ^/3i^mi' 



(1) 



with the order parameter ^(T*), defined on the subspace of the lowest Landau level (LLL)li3. This approximation 
is valid at sufficiently low temperatures or high magnetic fields ( i.e. for ksT <C hujc ), when thermal excitations 
to higher Landau levels in the condensate of Cooper pairs can be neglected. All possible configurations of the order 
parameter in this subspace can be taken into account by considering the free energy 

F ^-keT hi Z (2) 

with the partition function Z defined by the functional integral: 

Z= f D^pD^j* exp l-FGL/kBT] (3) 

In Eq.(||) a = a' {Hc2 — H) and /3 are phenomenological constants. The integral in Eq.(|^) is performed over all non 
equivalent states. 

An arbitrary wave function (in the symmetric gauge) from the LLL-subspace can be written as a one dimensional 
integral 



V'(x,y) = e"«y dgc(g)0,(x,2/) (4) 

where (f)q {x, y) = exp iqx ~ (y + f ) ^^ ^ Landau function with an orbital center located at —q/'i along the y-axis. 
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Note that all spatial lengths are measured in the units of the magnetic length an — [ch/eH] . A system of A''- 



vortices , with size L^ = clxvN along the x-direction , where ax is an arbitrary constant, is described by N coefficients 
c, = c (qi), q, = ^i with i = -N/2 + 1, ... , N/2. 

It is well known that the minimal value of the GL free energy functional is obtained when only ViV coefficients 
from the whole set of N coefficients are different from zero ( i.e. c^^/j^ 7^ , for n = —\fN/2 + 1, ..., \fN /2 ). At 
sufficiently low temperatures , when amplitude fluctuations are suppressed (see later) , this minimum corresponds to 
the minimum of the Abrikosov parameter: 

i5^=(l-^d\\i,\^/{^^d\\^.f^ (5) 

where V is the volume of the superconductor. From the definition of /3a it is seen that the absolute minimum, /?„ ~ I, 
is obtained for a spatially uniform order parameter. Any deviation from 1-01 = const leads to an increase in (3a- Under 
the constraint of the LLL subspace, however, ji/'l ^ can not be a constant ( since \ip\ = at the vortex cores ), and 
the minimum [3a — Pa — 1, 159 is obtained for the triangular Abrikosov lattice , which is the closest configuration to 
the homogeneous one. Other periodic lattices yield small positive deviations from j3a , while any departure of ^ from 
the quasi uniform distribution of the vortex lattice towards a localized structure leads to a drastic increase of the free 
energyEj. 

Thus, we conclude that the main correction to mean field order parameter arises from fluctuations of c„ and a^ , 
and take: 



^/N/2 

E 

-^/N/2+l 



C{q)= E -r.5{q-—) (6) 



a.i: 



Note that regardless of the choice of c„, tfj{x,y) is a periodic function of xwith a period a^- Therefore, the 
used representation of c (q) allows the order parameter to fluctuate only along the y-direction. Each coefficient, 
C71 — \cn\ e*'''" , describes a set of -y/iV Landau orbital centers , periodically arranged within a certain chain along the 
X -axis ('Bragg chain'). These 'Bragg chains' ( Fig. 1) are allowed to slide independently along their common axis, 
where the phase ipn, determines the relative position x^, = — (/?„/g„ of the n-th chain. The ideal lattice states are 
obtained by selecting c„ = cq exp (17^^) , < 7 < 7r/2c3. Then for an arbitrary rhombic lattice the lattice constant 
Qx in units of magnetic length and the angle 0/2 between the principal crystallographic axes are expressed through 

7 as a^ = tt/v 1 — {'^/tt) , cos0 = 7/7r. 

The partition function (0) can be therefore approximated by the multiple functional integral: 

^ ~ / n ^Cnrf< exp {-foL) , (7) 

where 

KsJ vAI — ^ — ' 

-^ n n,s,t 

with a = , "°^ ^ , B = —fE£^^-- \ — exp ( — ^ ). The inclusion of the factor vW in the denominator takes into 
account the fact that each term in the GL free energy functional , Eq.pl) (written in the discrete representation 
(jm) , corresponding to a certain Bragg chain in the vortex lattice, is degenerate ViV times. This degeneracy reflects 
the freezing of all ^/N internal degrees of freedom within a chain. As discussed above, by relaxing these degrees of 
freedom the free energy functional develops very large, highly unprobable fluctuations. Therefore, by dividing Fql 
with vTV" we are left with a single nondegenerate term , which represents the dominant degree of freedom for each 
chain. 

Note that the functional Eq. (|^) is invariant under an arbitrary linear shift of the phases (pn ■ 

ip'„ ^ (fn + an + b (9) 

This symmetry follows from the symmetry of the GL free energy Eq. (nh under the magnetic translation groupEll. 
The functional (^) has a set of local minima with a single amplitude for the entire set of lattice configurations: 

where 

sp 

is a reduced Abrikosov parameter , and z = 7r^/a^. These minima can be immediately obtained if we use the symmetry 
of the functional Jgl under the translation n ^ n + 1 . The Abrikosov parameter /3a , defined by (o) , is proportional 
to P^: (3a — \pKi3alax- The triangular lattice , with /3a = f3A (E3), corresponds to 7 = 7r/2, a^ = 2-Kl\fi. The mean 
field condensation energy per unit vortex is Fgl/N = — A", . 

Our main approximation at this point is based on the small value of the parameter A ~ e~'^, which enables us to 
neglect in Eq.(|^) all terms of the order higher than A^ , i.e. to retain , in addition to the first order terms in A, only 
the leading order terms in the phase ( (/?„) dependent part of the free energy. Thus, up to this order in A 

Igl = -S^|c„|^ + -^ l|c„|'* + 4A|c„|^|c„+i|^ + 4A^|c„_i| |c„+i| |c„|^cosxn (10) 

n n 

The angles Xn ^I'e linear combinations of the phases Lpn , 

X„ = -2(p„ + <y9„_l +(^„+l, (11) 

which are clearly invariant under the transformation (0). 



III. SHEAR MOTIONS : THE MELTING MECHANISM 

The model presented in the previous section is based on our observation that the low lying excitations of the 
Abrikosov vortex lattice are associated with the sliding motions of the lattice Bragg chains along, the principal 
crystallographic axes. These excitations are closely related to the soft shear modes discussed by Moorea in connection 
with the vortex lattice melting. i-. 

Let us consider this analogy in a greater detail. Following Ref.(u ) we first invoke a perturbative approach with 

respect to the Abrikosov vortex lattice solution c„ , by defining the displacements 6„ through c„ = c„ ' (l + bn) 
, and their normal mode (phonon) coordinates uj = j^t/j "^ {bn ^ b^) e^'^"' , with k = 4=^, I = — \/iV/2 \/N /2 , 

n 

and then expanding the free energy functional ( |l^) to second order in these coordinates. Omitting the details of 
calculations our result can be written as: 

^foL^^J2[p>^Kf + Q^h\'] (12) 



P0a 



fe>0 



where in the long wavelength limit k ^ Q , Pk ^ 2, and Qk -^ — A^fc*cos27. Thus, it can be readily shown that the 

relative variance of the order parameter diverges, i.e. : ((?/'^)~(^) ) / ("0) ^/fr- Therefore, we can conclude that 

although the lattice states with 7r/4 < 7 < tt/2 are thermodynamically stable {S/gl > 0), the soft modes Qk , h-.-^ , 
are responsible for infinite fluctuations of the order parameter. This result is similar to that obtained by MooreQ for a 
2D system of fluctuations. The divergence, which arises because of the perturbative nature of the above calculation, is 
stronger in our ID model. However, within the nonperturbative method developed below this divergence is removed. 
It should be noted that the soft mode described above is associated with the long wavelength component of the 
phase fluctuations in Eq.dTO); this can be seen by neglecting amplitude fluctuations, deflning fluctuating phases: 

(fif (n) = ipn- </?i^^ = ipn - -^T^n'^ 



and taking the continuous limit (see Eq.(ll)), i.e. Xn ^ ^ + d^f/dy , so that the relevant part in the free energy 
functional (^ can be written as: 

Sfph ^Ka f cos {n + d\f/dy^) dy « ]^Ka j {d\f/dy^f dy (13) 



where Ka ~ A^ ^ . It is instructive to compare this expression with that derived in (□) for the effective Hamiltonian 

associated with a smoothly varying phase {x, y), namely: Hph = hcQQ J d^r CS/^9) , where cge is an isotropic shear 
modulus of the vortex lattice, which is given approximately by 2%- 

The agreement between the two approaches is , however- incomplete , not only because of the one dimensional 
nature of our model (in contrast to the 2D analysis of Ref(cl)), but also because of the signiflcant difference in the 
'stiffness' parameters Ka, and cee , namely Ka/cgq ra X^ r^ 10^^. The reason for the disagreement can be understood 
within our approach by considering shear motions along families of Bragg chains with Miller indices higher than of 
the principal ones. For these families the values of Ux are relatively large , and the corresponding values of A^ are not 
small compared to unity. In the limit of very large Miller indices, a^; ^ 00 , and A^ ^ 1, so that the corresponding 
stiffness parameter approaches Cge, and becomes independeiit of the chain orientation, as in Moore's theory. 

Thus, in contrast to the isotropic shear model used in (□) , the appearance of the small parameter A^ in front of 
the leading phase dependent terms of the free energy functional ( |lO| ) implies that the shear motions along the two 
principal crystallographic axes cost a small fraction of the condensation energy , and so lead to signiflcant distortions 
of the vortex lattice along these particular directions at very low temperatures ( i.e. with respect to the mean held 

^^)- . . . 

Considering this low temperatures regime, the calculation of the partition function Z can be simplified considerably 

since amplitude fiuctuations can be neglected. The functional integrals in Eq.(g|) over the order parameter {ip,ip*} 

should be replaced by integrals over the new variables {c„, c*} = {|c„| , ipn}- 

For an arbitrarily large system due to the invariance of free energy under the magnetic translation group, the 

transformation Eq.(pT|) is degenerate, and the inverse transformation, ipn = <Pn {{Xm}), is not unique. It is determined 

up to a linear function of n. To determine the phases (/?„ uniquely , we have to impose additional conditions, 

determining the arbitrary constant a and b in Eq.(P). These conditions are equivalent to boundary conditions of 

Eq.(ni|). Linear boundary conditions lead to linear dependence of a and b on Xm- Since the determinant of any such 



transformation does not depend on the variables, the partition functions for various a and b differ by a constant factor. 
Therefore, instead of {(/5n} one can integrate over {xn} with the same free energy functional. The GL free energy has 
different values when the phases {xn} he within the interval [0,7r]. To satisfy this condition and to exclude double 
counting of fluctuation we integrate over interval Xn S [Oi'"'] ^nd allow the phases (p„ in Eq.(|ll|) to have arbitrary 
values. 

Omitting unimportant constant factor we obtain after integration over angle variables that the partition function 
can be written as: 

oo 

Z^zf^o, fll\cn\d\cn\e-f% (14) 

" 

where 

/, = ^j-a|c„|V|(|c„|V4A|c„|'|c„+i|') -iln/o(2/3A27r|c„_i||c„+i||c„|')| (15) 

and Ik (x) is the modified Bessel function of the order k. 

Neglecting amplitude fluctuations, the integrals in Eq.dlJ) can be performed by the stationary phase approximation. 
Since the last term in Eq.(|l5|) is of the order A^ or smaller, the approximate solution to the stationary point equations, 
Q l I = 0, can be simply obtained by using the translational symmetry of the free energy functional fs- It is similar 

to the mean field solution, |c„|" = ==— with the generalized Abrikosov parameter 

ppfi 

■^ ^l + 4X-i\^^ (16) 

where 

4A2 na^ Tcm 



(1 + 4A)' 2/3 T 

The temperature Tcm (cix) determines a smooth crossover from the mean field lattice state with 7 = 7r/2, (3 fi = 
/3; = 1 + 4A — 4A^, to a new state corresponding to Pn — /3„^ = 1 + 4A, where the phase dependent terms in the free 
energy are completely destroyed by fluctuations. Note that the energy difference between these states is of the order 
of the small parameter A^ . 

In the zero temperature limit, T ^ Tcm (a^), the parameter /3/; ~ ^ (l + 4A — 4A^) has minimal values at a^ = ^ 

, and a^, — 2^/3tt (Fig. 2) , depending on the choice of the 'Bragg chains' family (i.e. along the 2: or x axis in Fig. 
1). Both of the minima describe a triangular Abrikosov lattice with /3/; = Pa — 1.1596. Both directions can be 
selected in three equivalent ways in the Abrikosov lattice. All equivalent configurations can be obtained from the 
invariance of the mean field Abrikosov parameter f3a — \f^'Yli ^■'^P ("-^ (•^^ "^ P^)) '^'-'^ (^Tsp), where z = 7r^/a^, under 



the transformations z' — Ji/'-i , 7' = Jl^ 'i. and z' = 2;, 7' = —7 or 7' = 7 + Trn with an arbitrary integer n. 

The doubly degenerate , equilibrium state at T = 0, just described, is stabilized by the competition between two 
types of interactions among parallel chains: The repulsive interaction between any two neighboring chains, which is 
linear in the coupling parameter A , and the attractive 3-body, phase dependent interaction ( i.e. involving any three 



neighboring chains) , which is quadratic in A (see Eq.(lO)). At finite, low temperatures, i.e. when T ~ T^m , the 



shear fluctuations destroy the phase coherence among parallel Bragg chains, thus diminishing the small attractive 
interaction , and raising the total free energy. The relatively large, repulsive interaction is affected only at higher 
temperatures. 

2/2 

The interchain coupling parameter A , depends on the lattice parameter Ux, through A ~ e^'^ '"^ . Since a^' > ax 
(Fig. 1), the chains along x' are closer to each other than those along x , and so A {ax') > A (a^). Consequently, at 
low temperatures , T < Tcm (ax) , when the attractive 3-body interaction diminishes with increasing temperature, 
the flrst state {a^> ) is more stable than the second one (0:^), since its free energy increases more slowly with increasing 
temperature than that of the second one (Fig. 3). At higher temperatures T > Tcm (cix'), when the repulsive interchain 
couplings determine the temperature dependence, the tendency is reversed and the free energy of the flrst state {a^i ) 
increases faster with increasing temperature than that of the second one {ax)- Thus, there is an intersection point 
Tcm {ix) < Tm < Tcm {o-x') , at which the free energies of the these states are equal, but the corresponding entropies 



are a little different. Therefore, we conclude that at T = T,„ , there is a weak first order transition characterized 
by a small jump of the lattice entropy. Defining the parameter t = — g^"rp (EiJ), the position of the crossing point 
corresponds to t — i,„ ~ —16 , and the jump in the entropy ( S_= — T|y ) is AS" ~ 7.5 10~^Fmf/T . The values of 
tm and AS* agree pretty well with the Monte-Carlo simulationsM. 

The physical nature of this transition can be illuminated by considering the shear modulus /i. The vanishing of the 
shear modulus in atomic crystals is usually regarded as a definition of the crystal melting point. In our case fj, can be 
calculated by transforming c^^ = e*''" c„tJ and taking the limit: 

Note that the considered transformation shifts the phases Xn by —277. Therefore, the shear modulus is proportional 
to the phase factor in the free energy, i.e. /i oc (cosx„). Normalized by the mean field value, ^mf, where cosx„ = —1, 
it is reduced to 

M _h{r) (^^^ 



The dependence of the shear modulus on the parameter t is plotted in Fig. (4). At the transition point, t — tm, 
the value of the parameter a^ , corresponding to the minimum free energy, changes abruptly and the shear modulus 
jumps from /ii to fi2- It should be stressed that ^2 7^ 0. The residual shear energy on the high temperature side of 
the transition point reflects an incomplete melting at t = t^ ■ The 'liquid' state on this side of the transition point 
retains some degree of phase coherence between different chains, which continues to decrease gradually to zero with 
increasing temperature, reaching the complete liquid state only asymptotically. This behavior seems to be due to the 
persistence of long range periodic order along the chains axis in our model at any temperature. 

Interesting structural information on the 'quasi liquid' states described above can be obtained from the calculation 
of the average values 

(Xn) = —r-TT / c?X«X^exp(-TC0SXn) 
I'-'o (T) Ja 

with A: = 1, 2 . In the low and high temperature limits 

1/2 



(Xn)=^-V- {xl)=^'-(-) +- forr»l (18) 

iXn) = ^/2 + 2t/tt {xD - 7rV3 + 2t for r « 1 



The square root of the relative variance a = y (Xn) ~ iXn) / (Xn), is found to be cr ~ ^^ <C 1 in the low temperatures 

regime and a ~ l/-\/3 in the high temperatures one. These results show that with the temperature increase the 
fluctuations destroy the phase correlation between chains so that the SC state transforms from a frozen Abrikosov 
lattice at zero temperature to a new, 'liquid' state with strong vortex fluctuations. However, in contrast to the usual 
liquid state , here we flnd that the average vortex positions form a regular lattice with 7 — 7r/4 and tt^/o^ ~ 2.97, where 
a^ corresponds to the minimum free energy. For this lattice the angle, 0/2 , between the principal crystallographic 
axes corresponds to ~ 75°. 

The first order 'melting' point at t ^ tm thus corresponds to a discontinuous (rotational) symmetry breaking in 
this lattice of average vortex positions, from ~ 60° (Abrikosov lattice) on the low temperature side , to O ~ 75° 
on the high temperature one. 

IV. BRAGG-CHAINS PINNING AND THE ABSENCE OF LONG RANGE ORDER 

An intriguing issue in the theory of the vortex lattice melting concerns the existence or the absence of long range 
phase coherence in the SC mixed state. In this section we address the problem of long range order (LRO) and the 
related topic of vortex lines pinning , as they appear in our model. 

As discussed in Sec. (Ill), the fluctuating phases (/?„ can not be uniquely determined from Xnl they depend on the 
choice of the boundary conditions for Eq.(0). Since the general solution of the homogeneous equation (xn = 0) is 
a linear function of n (Eq. (0) , two constants (a and b) are required to uniquely determine (pn ■ A possible choice is 



to take a = b = , which corresponds to the selection (po = , and (p-i — ipi. The physical meaning of the first 
condition is that the chain , labeled n = ( i.e. located vertically at y = 0) , is pinned to a fixed 'horizontal' (i.e. 
along the x-axis) position. The second condition has a clear physical meaning in the long wavelength limit , namely 
that the horizontal displacement u^ = [86 /dy) of the vortex lines vanishes at the pinning site y = . 
The solution of Eq.(0) which satisfies this particular pinning condition is { for n > 0) : 



V±n =^{n-l) X±t + ^Xo 



1=0 



(19) 



This transformation enables us to calculate any correlation function of phase factors; in particular the pair corre- 
lation function: 

^0 Jo 



can be readily evaluated by using Eq.(p^) to yield: 



Io{-r) 



n+l 



(20) 



where the function v ^ v {v; n, n') — v {ly; n' , n) is defined ( for n > n') by : 2 1^ ~ ^'1 j 3't ly = ; \n — n'\ , for < 
V < n' ; |n — z^l . for n' < v < n ; and 0, for v > n. 

In the high temperature limit t ^ 1 , far above the melting point, the small argument expansion of the modified 
Bessel function I^, (— r) yields: 



(e'('^"'''^")) ex T^i^"-"'") ^ 5„,„, 



(21) 



meaning no phase correlation at all. 

In the low temperature limit r 3> 1, the asymptotic expansion of I^ (— t) leads ( for any n' > n ^ 1) to the 
expression: 



J^Vn.'-Vr^}\ 



pi7r(n' — n )/2 



exp 



2t 



{Anf 



(22) 



where An ^ n' — n, and n — n'/3 + 2n/3 — 1/2. Eq.(|22|) is identical to a second order cumulant expansion with 
i'Pn) = ^{Xn)n^ « 57m2. 

This result shows that a genuine long range phase correlation exists only at zero temperature; it also shows that a 
cluster of highly correlated chains can grow only along a pinned chain, since the phase fluctuations diverge with the 
distance from the pinning chain (see Eq.(IT9)) . Note that the position of such a chain is arbitrary since there is no 
energy cost to pinning in the GL theory used. In real samples the translational symmetry is broken by impurities , 
crystal defects , and the termination of the lattice at the sample surface, which can pin chains of orbital centers to 
fixed positions. A single pinning center located near a given chain may pin the entire chain due to the chain rigidity. 
To maximize the pinning strength, however, additional pining centers should be distributed uniaxially along the same 
chain, rather than randomly. 

Let us now study the range of SC order existing in_tJie vortex state at finite temperature. This can be done by 
considering the size dependence of the structure factorllil: 



^<^)-^ 



I{G 



where 



/(^) = fd^r\^j{r)\\' 



^■- 



b ^ — 



and G is a reciprocal lattice vector of the Abrikosov lattice with G^ — ^p^ , Gy = '^f^ — 2iybx 

bx = ixn) ax/2'K , and v , m integers. At zero temperature the long range order is reflected by the Bragg peaks 

with S [G] (X N . At finite temperature: 



Sit) 



'-e'"/' X 



(23) 



^ 5q^2^^Ci{n' + i',n,n',n + i')e '""-(" "'' 



9y 



where at r 3> 1 



s^ (p - s/3 + 1/3) 



exp 



2t 



+ i iXn) vin- n') 



v>0 



(24) 



with s = min {v, \n' — n\), p = max (i/, \n' — n\) ,ni = n' + v , n2 — n, n^ — n' , 77,4 = n + i^. Note that C4 is the 
four-chain phase correlation function appearing in the quartic term of the GL free energy. 

Now, since {xn) = X is independent of n , the sum over n yields a factor vTV , and Eg. (|23|) can be rewritten as: 



s(^\ ^Vn'- 



■'-e-^"l^ X 



\1\>V 



H(yx 



\1\<V 



byGy) 



exp 



iH^X 



byGy) 



l^l-oH/2^ 



/2r 



(25) 



This expression reflects the extreme anisotropy characterizing our Bragg chains model: Along the reciprocal lattice 



axis Gx ~ (i.e. for j/ = in Eq.(25) one finds perfect LRO, since 5 (Gx = 0, Gy) ^ N. For any Gx 7^ , however, 



the corresponding Bragg peaks reflect only the ID LRO within the real lattice chains, i.e. S {Gx ^ OjCy) ~ N^''^. 
This feature is due to the finite range of the off-diagonal phase correlation function ( i.e. to the absence of off-diagonal 
LRO ) at finite temperature. 

It should be stressed here that in the triangular Abrikosov lattice there are three equivalent ways to select the 
principal axes. Since the reciprocal lattice points with Gx = depend on our concrete choice of the coordinate system 
one may expect that by averaging over all three equivalent orientations, the size dependence of the structure factor 
will be isotropic , satisfying S ~ N'^ with 1/2 < ct < 1. |— ■ 

This result is consistent with the quasi LRO obtained by Kato et ahllil, which is reminiscent of the Kosterlitz- 
Thouless-Halperin-Nelson- Young (a) theory of 2D melting , according to which S ^ N°' with a < 5/6. 



V. CONCLUSION 



In this paper we have studied the melting of the SC vortex lattice in 2D at high magnetic fields and low temperatures 
(i.e. in the LLL approximation) by using a new approximate analytical approach. Our results basically agree with 
the state of the art Monte Carlo simulations. The simple analytical approach used enables us to draw a clear picture 
of the melting process: The skeleton of this picture consists of the principal crystallographic axes in the triangular 
Abrikosov lattice, along which families of almost rigid Bragg chains slide nearly freely at low temperatures due to 
thermal fiuctuations. Similar motions along crystal axes with higher Miller indices cost significantly more energy and 
are therefore quenched at low temperatures. 

The melting of the lattice occurs essentially when these fluctuations overcome the weak attractive interaction 
between chains. This interaction is not the same for the two principal axes; it is stronger for the more closely packed 
family of chains. 

Thus, the fluctuations destroy the order within the more loosely packed family of chains at lower temperature. 
Consequently the configuration based on the more closely packed family is the more stable one at low temperatures 
(i.e. below T„j ). However, at higher temperatures ( i.e. above Tm ) when the order in the closely packed family is 
also destroyed, the configuration based on the loosely packed family of chains becomes the more stable one , since the 
remaining interaction between chains at these temperatures is repulsive and weaker for the loosely packed family. 

The first order transition at Tm is therefore a discontinuous transformation between two different configurations 
of chains. The low temperature configuration, which is characterized by small fluctuations about the mean positions 
of vortices, forming an ideal Abrikosov triangular lattice (i.e. with angle 0/2 = 30° between the principal axes). In 
the high temperature configuration the phase fluctuations are significantly larger than in the low temperature one, 
whereas the mean positions of the vortices are still forming an exactly regular lattice, but now with « 75°. 



A correlated cluster of chains nucleates only around a pinned chain. Since according to our model, the pinning 
force of a whole chain can be strengthen dramatically by distributing pinning centers uniaxially along this chain, it 
may suggest a very efficient way of generating pitujiing defects in quasi 2D SC. This pinning mechanism may be tested 
experimentally by producing columnar defects (c3) along the conducting planes in quasi 2D SC. 
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Figure captions: 



a 



Fig.l The two families of Bragg chains in the triangular lattice along the principal axes x and x' . The parameters 
and Gx' are the periods of the order parameter modulus along these axes respectively, while ir/ax and it /ax' are 
the respective distances between chains. 

Fig. 2 Dependence of the Abrikosov parameter, /J^, on z = {tt/qx) ■ The two minima at zi = tt/2^/3 and Z2 = v37r/2 
correspond to the triangular Abrikosov lattice, (3a = Pa, with different choice of the Bragg chains direction. The 
maximum at z = tt/2 corresponds to the square lattice. 

Fig. 3 Free energy of fluctuating Bragg chains (normalized by the mean field free energy), —Pa/P/i, as a function of 
the reduced temperature, t = — "IZ^t ■ Solid and dash lines correspond to Bragg chains along the x and x' directions 
respectively. The intersection point at t ^ tm « — f6 determines the phase transition. 

Fig. 4 Dependence of the shear modulus, fi, on the reduced temperature, t. The jump oi ii at t — t„i reflects the 
'melting' transition. 
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